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In a gravitationally unstable universe, the structure of dark matter and galaxies, intergalactic 

^ , gas and magnetic field can have severe impact on the propagation of ultra high energy cosmic 

g ^ ■ rays (UHECR)Q The possible effects include spatial confinement and directional focusing along 

f-f ' the supergalactic matter sheets, as well as universal re-acceleration at large scale shock fronts, 

^^ ' and spectral modification due to energy dependent leakage into cosmic voids. As a result, the 

1/^ ' GZK-cutoff may be less pronounced and occur at a higher energy, where the stochastic nature 

of both acceleration and energy loss processes has to be taken into account. 
T-H . 
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^ ! 1 Supergalactic magnetic field structure and cosmic ray confinement 

("^ , Very little is known about the strength and orientation of the magnetic field outside 

^S| ' our galaxy. For cosmic ray transport calculations, one mostly uses the assumption 

^D , of a nanogauss field, which is homogeneous over cells with some reversal scale of 

^^ ■ order 1 Mpc. In such fields, the highest energy cosmic ray protons have a gyro-radius 

Tg ~ 300 Mpc, thus propagate almost in straight lines; this opens the possibility of 
an "UHECR astronomy", as anticipated by the Pierre Auger Projectul 

Models of structure formation in cosmology, however, draw a quite different pic- 
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Q ' ture: The magnetic field is aligned with the matter sheets, where it can reach a field 

|H , strength up to '^i/jG, while in the large cosmic voids the field drops to its primordial 

^ ' value of < pG;a this scenario is fully consistent with existing observations^ In the 

sheets, which have a typical thickness of ^lOMpc, the highest energy cosmic rays 
have Tg ^ 1 Mpc, and are thus confined. Outside the sheets, the accretion flow of in- 
k> , tergalactic gas drives the cosmic rays back, but the rapidly decreasing magnetic field 

;^ ■ may allow diffusive losses in upstream direction, which can imply spectral modifi- 

er 1 cations due to a "leaky box" mechanism. Fringe field effects may additionally focus 

and aUgn the cosmic rays with the field direction in the sheets; this might explain 
the apparent coocelation of UHECR arrival directions with the local sheet, the "super- 
galactic plane"QQ Since the universe needs no longer to be homogeneously filled with 
cosmic rays, the total energy budget for UHECR sources is strongly diminished. 

2 Large scale shocks and universal acceleration 

Another prediction in a structured universe is the existence of large scale shock frontsp 
providing the possibilityj-of cosmic ray acceleration by the very effective shock-drift 
acceleration mechanism^ In a global picture, the matter sheets form the collective 
downstream region, and the voids the collective upstream region in a foam-like shock 



topology. The cosmic rays, sliding sideways along the shock, never effectively leave 
the acceleration region. The spatial extension of the acceleration region in the di- 
rection of the flow can be estimated by the diffusion length, ^d, which depends on 
particle energy for quasi-perpendicular shocks and a Kolmogoroff turbulence spec- 
trum as lo ex E^/^. At the highest energies, ^d can become comparable to the sheet 
thickness, and the particles scatter freely between the boundary shock fronts. In this 
case, a stationary particle spectrum will no longer be obtained by the balance of dif- 
fusion over the shock front and downstream advection, but rather by the balance of 
energy gains and losses due to MBR interactions; here, the stochastic nature of the 
loss process turns out to be important. 

3 The stochastic nature of MBR pion production losses 

The transport of a proton in the MBR which is subject to pion production losses has 
to be described by a Markov point process, where the energy loss occurs randomly 
in distinct steps of random-distributed width. We may simplify the process to a pure 
counting process of unit steps, which is in case of a constant interaction rate known 
as a Poisson process. In photopion interactions, particles lose energy fractionally, i.e 
AE/E « A\nE — const « 0.2. For a Poisson process, one can show that an initial 
spectrum power law spectrum, / ex i?^°, of a source at distance D, suffers an energy 
independent reduction by a factor M — exp[— (D/A)(l — e^")], if A is the mean 
interaction length and a = a{A\nE). For a linearly increasing interaction rate, p ~ 
c/X = p'ln{E/Eo), the modified spectrum can_be approximately described as a 
power law steepened by Aa — [Dp' /c){l — e^")Ll 

The interaction rate in the microwave background can be best modeled relative 
to the maximum rate pi, which is reached for E > Ei^l ZeV and corresponds 
to Ai «4Mpc. For i?o«30EcV < E < Ei it is Hnearly increasing, (p/pi) ~ 
0.3\n{E/Eo), and p = for E < Eq. A continuing initial power law spectrum 
maps then to a piecewise power law with index a for E < Eq and E > Ei, and 
a + Aa in between. A spectral cutoff in the source maps to an exponential decline 
of the observed spectrum somewhat below the source cutoff energylZl We may give 
two numerical examples: A radio galaxy at D — 30Mpc, producing a spectrum 
/ oc E^^ with a sharp cutoff at 1 ZeV, is observed with a power law index a'~2.75 
between 30 and 300 EeV, followed by an exponential cutoff. A topological defect at 
D — 100 Mpc, producing a / oc E~^-^ spectrum, is observed with a power law index 
a' ^3 between 30 EeV and 1 ZeV, flattening back to a=2 for higher energies. 

4 Consequences for the GZK cutoff and cosmic ray observatories 

The time scale of large scale shock acceleration, ta, is generally larger than the time 
scale for MBR photopion losses, t^^; depending on magnetic field strength and shock 



velocity, we may find ratios ia/^Tr ~ 1 — 100 at >100 EeVH Thus the acceleration is 
not really effective in the ordinary sense; however, considering the stochastic nature 
of energy losses and the breakdown of advection at the highest energies, the resulting 
stationary spectrum can still be relatively flat for ii^>100 EeV: In the simple case of a 
constant acceleration time scale, interaction loss balanced shock acceleration leads to 
power laws / oc E^^, and the relation ta/iir = ^'[1 — exp(— 0.26)]^^ holds. Spectral 
indices as observed in the UHECR spectrum are obtained for ta/iyr ~ 4, but steepen 
very fast for larger values; under realistic conditions, the equilibrium spectrum is 
probably concave and too steep to explain the highest energy event rates. 

Therefore, the existence of large scale shocks in the universe does not make cos- 
mic ray point sources unnecessary; radio galaxies, AGN, gamma ray bursters or topo- 
logical defects may still contribute as UHECR sources. Clusters of galaxies, which 
are the sites of the strongest large scale shocks and well located in the universe, can 
play an intermediate role between point sources and large scale acceleration^ The im- 
portance of large scale shocks is rather that they provide a re-acceleration mechanism 
which is as universal as the GZK process, and thus may lead to revised estimates of 
the maximum distance of the possible sources of highest energy cosmic rays. Conse- 
quently, the pros and cons for the various source models have to be reconsidered in 
a structured universe. For the Pierre Auger UHECR observatory, the large values of 
the magnetic field arising from large scale structure simulations give little hope to see 
point sources of charged cosmic ray particles; however, UHECR events are likely to 
occur in clusters and map the local large scale structure of the magnetic field. 
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